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Abstract 

This paper investigates the problem of determining a binary-valued function through a sequence of strategically 
selected queries. The focus is an algorithm called Generalized Binary Search (GBS). GBS is a well-known greedy 
algorithm for determining a binary-valued function through a sequence of strategically selected queries. At each step, 
a query is selected that most evenly splits the hypotheses under consideration into two disjoint subsets, a natural 
generalization of the idea underlying classic binary search. This paper develops novel incoherence and geometric 
conditions under which GBS achieves the information-theoretically optimal query complexity; i.e., given a collection 
of TV hypotheses, GBS terminates with the correct function after no more than a constant times log TV queries. 
Furthermore, a noise-tolerant version of GBS is developed that also achieves the optimal query complexity. These 
results are applied to learning halfspaces, a problem arising routinely in image processing and machine learning. 



I. Introduction 

This paper studies learning problems of the following form. Consider a finite, but potentially very large, collection 
of binary-valued functions T-L defined on a domain X. In this paper, T-L will be called the hypothesis space and X 
will be called the query space. Each h E H is a mapping from X to { — 1, 1}. Throughout the paper we will let TV 
denote the cardinality of H. Assume that the functions in H are unique and that one function, h* E H, produces 
the correct binary labeling. It is assumed that h* is fixed but unknown, and the goal is to determine h* through as 
few queries from X as possible. For each query x E X, the value h*{x), possibly corrupted with independently 
distributed binary noise, is observed. The goal is to strategically select queries in a sequential fashion in order to 
identify h* as quickly as possible. 

If the responses to queries are noiseless, then the problem is related to the construction of a binary decision tree. 
A sequence of queries defines a path from the root of the tree (corresponding to T-L) to a leaf (corresponding to a 
single element of %). There are several ways in which one might define the notion of an optimal tree; e.g., the tree 
with the minimum average or worst case depth. In general the determination of the optimal tree (in either sense 



above) is a combinatorial problem and was shown by Hyafil and Rivest to be NP-complete |19|. Therefore, this 
paper investigates the performance of a greedy procedure called generalized binary search (GBS), depicted below 
in Fig. [T] At each step GBS selects a query that results in the most even split of the hypotheses under consideration 
into two subsets responding +1 and —1, respectively, to the query. The correct response to the query eliminates 
one of these two subsets from further consideration. We denote the number of hypotheses remaining at step n by 
\Hn\- The main results of the paper characterize the worst-case number of queries required by GBS in order to 
identify the correct hypothesis h* . More formally, we define the notion of query complexity as follows. 
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Generalized Binary Search (GBS) 

initialize: n = 0, Ho = H. 
while \Hn\>l 

1) Select Xn = argmiiia^gAr | J2heHn 

2) Query with Xn to obtain response t/„ = h*{xn)- 

3) Set Un+l = {h&Un'- h{Xn) = Un}, U = U + 1. 

Fig. 1. Generalized Binary Searcli, also known as the Splitting Algorithm. 

Definition 1: The minimum number of queries required by GBS (or another algorithm) to identify any hypothesis 
in "H is called the query complexity of the algorithm. The query complexity is said to be near-optimal if it is within 
a constant factor of logN, since at least logN queries are required to specify one of N hypotheses. 

Conditions are estabUshed under which GBS (and a noise-tolerant variant) have a near-optimal query complexity. 
The main contributions of this paper are two-fold. First, incoherence and geometric relations between the pair 
(Af, 'H) are studied to bound the number of queries required by GBS. This leads to an easily verifiable sufficient 
condition that guarantees that GBS temiinates with the correct hypothesis after no more than a constant times 
logN queries. Second, noise-tolerant versions of GBS are proposed. The following noise model is considered. 
The binary response y G {—1,1} to a query x G X is m independent realization of the random variable Y 
satisfying F{Y = h*{x)) > P(y = —h*{x)), where P denotes the underlying probability measure. In other words, 
the response to x is only probably correct. If a query x is repeated more than once, then each response is an 
independent realization of Y. A new algorithm based on a weighted (soft-decision) GBS procedure is shown to 
confidently identify h* after a constant times logN queries even in the presence of noise (under the sufficient 
condition mentioned above). An agnostic algorithm that performs well even if h* is not in the hypothesis space % 
is also proposed. 

A. Notation 

The following notation will be used throughout the paper. The hypothesis space 'H is a finite collection of binary- 
valued functions defined on a domain X, which is called the query space. Each /i e 'H is a mapping from X to 
{—1,1}. For any subset %' d H, \H'\ denotes the number of hypotheses in %' . The number of hypotheses in % 
is denoted by N := |^|. 

II. A Geometrical View of Generalized Binary Search 

The efficiency of classic binary search is due to the fact at each step there exists a query that splits the pool of 
viable hypotheses in half. The existence of such queries is a result of the special ordered structure of the problem. 
Because of ordering, optimal query locations are easily identified by bisection. In the general setting in which the 
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query and hypothesis space are arbitrary it is impossible to order the hypotheses in a similar fashion and "bisecting" 
queries may not exist. For example, consider hypotheses associated with halfspaces of X = Each hypothesis 
takes the value +1 on its halfspace and —1 on the complement. A bisecting query may not exist in this case. To 
address such situations we next introduce a more general framework that does not require an ordered structure. 



A. Partition of X 

While it may not be possible to naturally order the hypotheses within X, there does exist a similar local geometry 
that can be exploited in the search process. Observe that the query space X can be partitioned into equivalence 
subsets such that every h E H is constant for all queries in each such subset. Let A{X,H) denote the smallest 
such partitiorQ Note that X = UAeyi^- ^'^^ every A £ A and h £ H, the value of h{x) is constant (either +1 
or —1) for all x £ A; denote this value by h{A). Observe that the query selection step in GBS is equivalent to an 
optimization over the partition cells in A. That is, it suffices to select a partition cell for the query according to 
An = argmiuAe^ | Y.heH„ 

The main results of this paper concern the query complexity of GBS, but before moving on let us comment on 
the computational complexity of the algorithm. The query selection step is the main computational burden in GBS. 
Constructing A may itself be computationally intensive. However, given A the computational complexity of GBS 
is N \A\, up to a constant factor, where \A\ denotes the number of partition cells in A. The size and construction 
of A is manageable in many practical situations. For example, if X is finite, then |^| < \X\, where \X\ is the 



cardinality of X. Later, in Section V-B we show that if % is defined by N halfspaces of X := M^, then \A\ grows 
like N'^. 



B. Distance in X 

The partition A provides a geometrical link between X and T-L. The hypotheses induce a distance function on A, 
and hence X. For every pair A,A'eA the Hamming distance between the response vectors {hi{A), . . . , h^iA)} 
and {hi{A'), . . . , hN{A')} provides a natural distance metric in X. 

Definition 2: Two sets A, A' E A are said to be k-neighbors if k or fewer hypotheses (along with their 
complements, if they belong to %) output different values on A and A' . 

For example, suppose that T-L is symmetric, so that h E H implies —heT-L. Then two sets A and A' are k- 
neighbors if the Hamming distance between their respective response vectors is less than or equal to 2k. If H is 
non-symmetric (h £ H implies that —h is not in H), then A and A' are fc-neighbors if the Hamming distance 
between their respective response vectors is less than or equal to k. 

'Each h splits X into two disjoint sets. Let := {x a X : h(x) = +1} and let Cn denote its complement. A is the collection of all 
non-empty intersections of the form nhgir '^h. where £ {Cjj, (5^}, and it is the smallest partition that refines the sets {ChjhgH- is 
known as the join of the sets {Ch^hsT-L- 
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Definition 3: The pair {X, H) is said to be k-neighborly if the fc-neighborhood graph of A is connected (i.e., 
for every pair of sets in A there exists a sequence of fc-neighbor sets that begins at one of the pair and ends with 
the other). 

If {X, %) is fc-neighborly, then the distance between A and A' is bounded by k times the minimum path length 
between A and A' . Moreover, the neighborly condition implies that there is an incremental way to move from one 
query to the another, moving a distance of at most k at each step. This local geometry guarantees that near-bisecting 
queries almost always exist, as shown in the following lenmia. 

Lemma 1: Assume that {X,'hL) is k-neighborly and define the coherence parameter 



c*(X,T-L) := minmax 
p hen 



J2h{A)P{A) , (1) 
AeA 

where the minimization is over all probability massfimctions on A For every H' C7i and any constant c satisfying 
c* < c < 1 there exists an A G A that approximately bisects H' 



hen' 



< c\H'\ 



or the set %' is a small 



\H'\ < 



k 



where \'H'\ denotes the cardinality ofH'. 

Proof: According to the definition of c* it follows that there exists a probability distribution P such that 



E E h{A)P{A) 
hew AeA 

This implies that there exists an A e ^ such that 



< c\U'\ 



hew 



< c\n'\ 



or there exists a pair A and A' such that 

E HA) > c\H'\ and ^ H^) < "cl^'l • 



hew 



hew 



In the former case, it follows that a query from A will reduce the size of %' by a factor of at least (1 + c)/2 (i.e., 
every query x & A approximately bisects the subset %'). In latter case, an approximately bisecting query does not 
exist, but the fc-neighborly condition implies that \'H'\ must be small. To see this note that the ^-neighborly condition 
guarantees that there exists a sequence of fc-neighbor sets beginning at A and ending at A' . By assumption in this 
case, > on every set and the sign of "^hew H') must change at some point in the sequence. 

It follows that there exist fc-neighbor sets A and A' such that J^hew > cl'H'l and '^hew H^^') < 

—c\H'\. Two inequaUties follow from this observation. First, J2heW H-^) ~ J2heW H-^') > ^c\H'\. Second, 
I X^/jg^/ ^(^) ~ T^heW '^(^01 — Note that if h and its complement h' belong to %' , then their contributions 
to the quantity | J^hew ~ J2hen' H-^')] cancel each other. Combining these inequalities yields \H'\ < k/c. 
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Example 1: To illustrate Lemma^ consider the special situation in which we are given two points Xi,X2 G X 
known to satisfy h*{xi) = +1 and h* {X2) = —1. This allows us to restrict our attention to only those hypotheses 
that agree with h* at these points. Let % denote this collection of hypotheses. A depiction of this situation is shown 
in Fig. [2] where the solid curves represent the classification boundaries of the hypotheses, and each cell in the 
partition shown corresponds to a subset of X (i.e., an element of A). As long as each subset is non-empty, then 
the 1-neighborhood graph is connected in this example. The minimization in (|7]l is achieved by the distribution 
P = \^xi + \6x2 (equal point-masses on xi and X2) and c*{X,'H) — 0. Lemma^implies that there exists a 
query (equivalently a partition cell A) where half of the hypotheses take the value +1 and the other half —1. The 
shaded cell in Fig. [2] has this bisection property. The figure also shows a dashed path between xi and X2 that 
passes through the bisecting cell. 




Fig. 2. An illustration of the idea of GBS. Each solid curve denotes the decision boundary of a hypothesis. There are six boundaries/hypotheses 
in this example. The connect hypothesis in this case is known to satisfy h*(xi) = +1 and h*(x2) = — 1. Without loss of generality we may 
assume that all hypotheses agree with h* at these two points. The dashed path between the points xi and X2 reveals a bisecting query location. 
As the path crosses a decision boundary the con'esponding hypothesis changes its output from +1 to —1 (or vice- versa, depending on the 
direction followed). At a certain point, indicated by the shaded cell, half of the hypotheses output +1 and half output —1. Selecting a query 
from this cell will bisect the collection of hypotheses. 

C. Coherence and Query Complexity 

The coherence parameter c* quantifies the informativeness of queries. The coherence parameter is optimized over 
the choice of P, rather than sampled at random according to a specific distribution on X, because the queries may 
be selected as needed from X. The minimizer in ([T]i exists because the minimization can be computed over the 
space of finite-dimensional probability mass functions over the elements of A. For c* to be close to 0, there must 
exist a distribution P on ^ so that the moment of every h E H is close to zero (i.e., for each h E H the probabilities 
of the responses +1 and —1 are both close to 1/2). This implies that there is a way to randomly sample queries 
so that the expected response of every hypothesis is close to zero. In this sense, the queries are incoherent with 
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the hypotheses. In Lemma [T] c* bounds the proportion of the split of any subset %' generated by the best query 
(i.e., the degree to which the best query bisects any subset %'). The coherence parameter c* leads to a bound on 
the number of queries required by GBS. 

Theorem 1: If [X ^T-L) is k-neighborly, then GBS terminates with the correct hypothesis after at most 
[log7V/log(A~^)] queries, where A = max{i±^, |^}- 

Proof: Consider the rtth step of the GBS algorithm. Lemma [T] shows that for any c e [c*, 1) either there exists 
an approximately bisecting query and \'Hn\ < or l^n-il < fc/c. The uniqueness of the hypotheses 

with respect to X implies that there exists a query that eliminates at least one hypothesis. Therefore, < 
— 1 = ["Hji-iKl — I'Hji-il^^) < ["Hn-iKl — c/k). It follows that each GBS query reduces the number of 
viable hypotheses by a factor of at least 

A :— minmaxl — ^, 1 — c/fc 1 = max 
c>c* 12 J 

Therefore, \Hn\ < iVA" and GBS is guaranteed to terminate when n satisfies A^A" < 1. Taking the logarithm of 
this inequality produces the query complexity bound. ■ 
Theorem [l] demonstrates that if {X ,%) is neighborly, then the query complexity of GBS is near-optimal; i.e., 
within a constant factor of log2 N . The constant depends on coherence parameter c* and fc, and clearly it is desirable 
that both are as small as possible. Note that GBS does not require knowledge of c* or k. We also remark that the 
constant in the bound is not necessarily the best that one can obtain. The proof involves selecting c to balance 
splitting factor and the "tail" behavior 1 — c/k, and this may not give the best bound. The coherence parameter 
c* can be computed or bounded for many pairs {X, Ti) that are commonly encountered in applications, as covered 
later in Section [V] 

III. Noisy Generalized Binary Search 

In noisy problems, the search must cope with erroneous responses. Specifically, assume that for any query x € X 
the binary response y e {—1, 1} is an independent realization of the random variable Y satisfying P(F = h*{x)) > 
¥{Y — ~h*{x)) (i.e., the response is only probably correct). If a query x is repeated more than once, then each 
response is an independent reaUzation of Y. Define the noise-level for the query x as '■— P{Y ~ —h*{x)). 
Throughout the paper we will let a :— sup^g;^- and assume that a < 1/2. Before presenting the main approach 
to noisy GBS, we first consider a simple strategy based on repetitive querying that will serve as a benchmark for 
comparison. 

A. Repetitive Querying 

We begin by describing a simple noise-tolerant version of GBS. The noise-tolerant algorithm is based on the 
simple idea of repeating each query of the GBS several times, in order to overcome the uncertainty introduced 
by the noise. Similar approaches are proposed in the work Kaariainen p6) . Karp and Kleinberg pO) analyze 
of this strategy for noise-tolerant classic binary search. This is essentially like using a simple repetition code to 



j l + c* k + l \ 
1 2 ' fc + 2j 
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communicate over a noisy channel. This procedure is termed noise-tolerant GBS (NGBS) and is summarized in 
Fig. [3] 

Noise-Tolerant GBS (NGBS) 

input: H, repetition rate R> 1 (integer), 
initialize: n = 0, Ho ~ H. 
while |-H„| > 1 

1) Select A„ = aTgrniuAeA I Y.hen„ 

2) Query R times from An to obtain R noisy versions 
of Un = h*{An). Let yn denote the majority vote of 
the noisy responses. 

3) Set Hn+i ={h&Hn: h{An) = y„}, n = n + 1. 

Fig. 3. Noise-tolerant GBS based on repeated queries. 



Theorem 2: Let hq denote the number of queries made by GBS to determine h* in the noiseless setting. Then in 
the noisy setting, with probability at least max{0 , 1 — rtoe^^'^""! } the noise-tolerant GBS algorithm in Fig. [i] 
terminates in exactly Rhq queries and outputs h*. 

Proof: Consider a specific query x G A" repeated R times, let p denote the frequency of +1 in the R trials, 
and let p = E[p]. The majority vote decision is correct if |p — p| < ^ — a. By Chernoff's bound we have 
IP(|p ^ ^ ~ Oi) < 2e^^^l3^"l'. The results follows by the union bound. ■ 

Based on the bound above, R must satisfy R > to guarantee that the labels determined for all tiq queries 

are correct with probability 1 — S. The query complexity of NGBS can thus be bounded by "ii'/^-ai^*^^ ■ Recall 
that N = \T-L\, the cardinality of H. If no = logiV, then bound on the query complexity of NGBS is proportional 
to log log a logarithmic factor worse than the query complexity in the noiseless setting. Moreover, if an 

upper bound on ng is not known in advance, then one must assume the worst-case value, rio = N, in order to set R. 
That is, in order to guarantee that the correct hypothesis is determined with probability at least 1 — 6, the required 
number of repetitions of each query is i? = [ |'°/2^q'|2 1 ■ In this situation, the bound on the query complexity of 
NGBS is proportional to log log ^, compared to log in the noiseless setting. It is conjectured that the extra 
logarithmic factor cannot be removed from the query complexity (i.e., it is unavoidable using repetitive queries). 
As we show next, these problems can be eliminated by a more sophisticated approach to noisy GBS. 

B. Soft-Decision Procedure 

A more effective approach to noisy GBS is based on the following soft-decision procedure. A similar procedure 
has been shown to be near-optimal for the noisy (classic) binary search problem by Burnashev and Zigangirov 
||3) and later independently by Karp and Kleinberg pO) . The crucial distinction here is that GBS calls for a 
more general approach to query selection and a fundamentally different convergence analysis. Let po be a known 
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probability measure over H. That is, pq : H ^ [0, 1] and J^henPoW ~ 1- The measure pq can be viewed as an 
initial weighting over the hypothesis class. For example, taking pq to be the uniform distribution over H expresses 
the fact that all hypothesis are equally reasonable prior to making queries. We will assume that po is uniform for 
the remainder of the paper, but the extension to other initial distributions is trivial. Note, however, that we still 
assume that h* E H is fixed but unknown. After each query and response {xn,yn), n — 0,1, ... , the distribution 
is updated according to 

p„+i(M oc p„(/i)/3(i--(''))/^(l (2) 

where Zn(h) = h{xn)yn, ft. G /3 is any constant satisfying < /? < 1/2, and pn+i{h) is normalized to 
satisfy '^^heuPn+iW — 1 ■ The update can be viewed as an application of Bayes rule and its effect is simple; 
the probability masses of hypotheses that agree with the label y„ are boosted relative to those that disagree. The 
parameter /3 controls the size of the boost. The hypothesis with the largest weight is selected at each step: 

hn := argm;«p„(/i). 

If the maximizer is not unique, one of the maximizers is selected at random. Note that, unlike the hard-decisions 
made by the GBS algorithm in Fig. [T] this procedure does not eliminate hypotheses that disagree with the observed 
labels, rather the weight assigned to each hypothesis is an indication of how successful its predictions have been. 
Thus, the procedure is termed Soft-Decision GBS (SGBS) and is summarized in Fig. [4] 

The goal of SGBS is to drive the error P(/i„ 7^ h*) to zero as quickly as possible by strategically selecting 
the queries. The query selection at each step of SGBS must be informative with respect to the distribution p„. 
In particular, if the weighted prediction '^ZheH Pri{h)h{x) is close to zero for a certain x (or A), then a label 
at that point is informative due to the large disagreement among the hypotheses. If multiple A E A minimize 
I X]/i6«P«(^)^(^)l' '^hen one of the minimizers is selected uniformly at random. 

Soft-Decision Generalized Binary Search (SGBS) 

initialize: pq uniform over %. 
forn==0,l,2,... 

1) An = argmin^e^ I I],jg^p„(ft)ft(A)|. 

2) Obtain noisy response j/„. 

3) Bayes update pn — > Pn+i\ Eqn. (|2]l. 
hypothesis selected at each step: 

K ■= argmax,ig//p„(/i) 



Fig. 4. Soft-Decision algoritlim for noise-tolerant GBS. 

To analyze SGBS, define C„ :— (1 — Pn{h*)) /pn{h*), n>0. The variable C„ was also used by Burnashev and 
Zigangirov |[3| to analyze classic binary search. It reflects the amount of mass that p„ places on incorrect hypotheses. 
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Let P denotes the underlying probability measure governing noises and possible randomization in query selection, 
and let E denote expectation with respect to P. Note that by Markov's inequality 

nK ^ h*) < P(p„(/i*) < 1/2) = P(C„ > 1) < E[C„] . (3) 

At this point, the method of analyzing SGBS departs from that of Burnashev and Zigangirov Q which focused only 
on the classic binary search problem. The lack of an ordered structure calls for a different attack on the problem, 
which is summarized in the following results and detailed in the Appendix. 

Lemma 2: Consider any sequence of queries {x„}„>o and the corresponding responses {yn}n>o- U P ^ ^> 
then {C„}n>o nonnegative supermartingale with respect to {pn}n>o; i-^-, ^[Cn+i\Pn] < Cn for all n > 0. 

The lemma is proved in the Appendix. The condition /? > a ensures that the update (j2]) is not overly aggressive. 
It follows that E[C„] < Co and by the Martingale Convergence Theorem we have that lim„^oo Cn exists and is 
finite (for more information on Martingale theory one can refer to the textbook by Bremaud |37|). Furthermore, 
we have the following theorem. 

Theorem 3: Consider any sequence of queries {a;„}„>o cind the corresponding responses {j/n}n>o- If P > a, 
then lim„_>.oo P(^ri ^ h* ) = Q. 

Proof: First observe that for every positive integer n 

E[C„] = n{Cn/C.n^i)Cn^i] = E[E[(C„/C„_i)C„-ib„-i]] 

= E[C„_iE[(C„/C„„i)|p„_i]] < E[C„-i] maxE[(C„/C„_i)b„_i] 

< Co ( max max E[(Ci+i/Ci)|pi] 

\ z— 0,...,n— 1 Pi 

In the proof of Lemma |2] it is shown that if /3 > a, then E[(Ci+i/Ci)|pi] < 1 for every pi and therefore 
maxp. E[(Ci+i/Ci)|pi] < 1. It follows that the sequence a„ :— (maxi^o....,n-i maxp. E[(Ci+i/Ci)|pi])" is 
monotonically decreasing. The result follows from ■ 

Note that if we can determine an upper bound for the sequence maxp^ E[(Ci+i/Ci)|pi] < 1 — A< 1, i = 
0, . . . ,n - 1, then it follows that that P(/J„ ^ h*) < N{1 - A)" < iVe^^". Unfortunately, the SGBS algorithm 
in Fig. |4] does not readily admit such a bound. To obtain a bound, the query selection criterion is randomized. 
A similar randomization technique has been successfully used by Burnashev and Zigangirov Q and Karp and 
Kleinberg | 30| to analyze the noisy (classic) binary search problem, but again the lack of an ordering in the general 
setting requires a different analysis of GBS. 

The modified SGBS algorithm is outlined in Fig. |5] It is easily verified that Lemma |2] and Theorem [3] also hold 
for the modified SGBS algorithm. This follows since the modified query selection step is identical to that of the 
original SGBS algorithm, unless there exist two neighboring sets with strongly bipolar weighted responses. In the 
latter case, a query is randomly selected from one of these two sets with equal probability. 

For every A E A and any probability measure p onH the weighted prediction on A is defined to be W{p, A) := 
X^/ieff where h{A) is the constant value of h for every x E A. The following lemma, which is the 

soft-decision analog of Lemma [T] plays a crucial role in the analysis of the modified SGBS algorithm. 
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Modified SGBS 

initialize: pq uniform over H. 
for n = 0,1,2,... 

1) Let b = min^g^ I ^^g^p„(/i)/i(v4)|. If there exist 1-neighbor sets A and A' with 
^^g^ p„(/i)/i(A) > b and J2he'HPn(^)^(^') ~^ ' '-'^^'^ select a;„ from A or A' with 
probability 1/2 each. Otherwise select Xn from the set A,,,,,, argmin^gyi | X]/ig-H^'"(^)^(^)l- ^'^ 
the case that the sets above are non-unique, choose at random any one satisfying the requirements. 

2) Obtain noisy response ?/„. 

3) Bayes update p„ — > Pn+i, Eqn. 
hypothesis selected at each step: 
hn := aigmaxheH Pnih) 

Fig. 5. Modified SGBS Algorithm. 



Lemma 3: If [X ^T-L) is k-neighborly, then for every probability measure p on % there either exists a set A Cz A 
such that \W{p, A)\ < c* or a pair of k-neighbor sets A, A' £ A such that W{p, A) > c* and W{p, A') < — c*. 

Proof: Suppose that min^g,^ A)| > c*. Then there must exist A, A' G A such that W{p,A) > c* 

and W{p,A') < — c*, otherwise c* cannot be the incoherence parameter of H, defined in ([T]l. To see this 
suppose, for instance, that W{p, A) > c* for all A E A. Then for every distribution P on A" we have 
J2AeA^he'HPi^)^i^)-^i^) > c*. This contradicts the definition of c* since J2AeA^he'HPi^)^i^)^i^) — 
X^/ie-H X^Ae^ '^'C^) — ™9'X'iG'H I X^Aeyt '*(^) -^(^)l- neighborly condition guarantees that there 

exists a sequence of fc-neighbor sets beginning at A and ending at A'. Since A)| > c* on every set and the 

sign of W{p,-) must change at some point in the sequence, it follows that there exist fc-neighbor sets satisfying 
the claim. ■ 

The lemma guarantees that there exists either a set in A on which the weighted hypotheses significantly disagree 
(provided c* is significantly below 1) or two neighboring sets in A on which the weighted predictions are strongly 
bipolar. In either case, if a query is drawn randomly from these sets, then the weighted predictions are highly variable 
or uncertain, with respect to p. This makes the resulting label informative in either case. If {X ,H) is 1-neighborly, 
then the modified SGBS algorithm guarantees that P(/i„ 7^ ft.*) — > exponentially fast. The 1-neighborly condition 
is required so that the expected boost to Pn{h*) is significant at each step. If this condition does not hold, then 
the boost could be arbitrarily small due to the effects of other hypotheses. Fortunately, as shown in Section [V] the 
1-neighborly condition holds in a wide range of common situations. 

Theorem 4: Let P denote the underlying probability measure (governing noises and algorithm randomization). If 
/3 > a and (X , H) is 1-neighborly, then the modified SGBS algorithm in Fig. pjgenerates a sequence of hypotheses 
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satisfying 

¥{K^h*) < N{l-Xy' < TVe-^" , n = 0, 1,... 

with exponential constant A = min | ^^"^ , || ^1 — — )' ^^''^''^ ^* '■^ defined in (1 \. 

The theorem is proved in the Appendix. The exponential convergence rat^is governed by the coherence parameter 
< c* < 1. As shown in Section [V], the value of c* is typically a small constant much less than 1 that is independent 
of the size of H. In such situations, the query complexity of modified SGBS is near-optimal. The query complexity 
of the modified SGBS algorithm can be derived as follows. Let (5 > be a pre-specified confidence parameter. The 
number of queries required to ensure that P(/i„ ^ h*) < S is n > X^^ log ^, which is near-optimal. Intuitively, 
about logiV bits are required to encode each hypothesis. More formally, the noisy classic binary search problem 
satisfies the assumptions of Theorem |4] (as shown in Section V-A i, and hence it is a special case of the general 



problem. Using information-theoretic methods, it has been shown by Burnashev and Zigangirov (3) (also see the 
work of Karp and Kleinberg 130)) that the query complexity for noisy classic binary search is also within a constant 



factor of log ^ . In contrast, the query complexity bound for NOBS, based on repeating queries, is at least logarithmic 
factor worse. We conclude this section with an example applying Theorem |4] to the halfspace learning problem. 
Example 2: Consider learning multidimensional halfspaces. Let X = and consider hypotheses of the form 

hi{x) := sign{{a,, x) + bi) . (4) 

where Ui G M.'^ and foj G M parameterize the hypothesis hi and {a.i , x) is the inner product in M'*. The following 
corollary characterizes the query complexity for this problem. 

Corollary 1: Let % be a finite collection of hypotheses of form Q and assume that the responses to each query 
are noisy, with noise bound a < 1/2. Then the hypotheses selected by modified SGBS with (3 > a satisfy 

., > 1 B(l-a) a(l~B)\ 

Wl 

Section \ V-B\ The polynomial-time computational complexity follows from the work of Buck ^42| /, as discussed in 
Suppose that 7i is an e-dense set with respect to a uniform probability measure on a ball in M.'^ (i.e., 



th \ = ^ [l — ^''l^p^ ^ "^^ff ^'' j" Moreover, hn can be computed in time polynomial in \H\. 

The error bound follows immediately from Theorem^since c* — and is 1-neighborly, as shown in 



Section 



V-B 



for any hyperplane of the form Q TL contains a hypothesis whose probability of error is within e of it). The size 
of such an % satisfies logN < Cdloge"^, for a constant C > 0, which is the proportional to the minimum 
query complexity possible in this setting, as shown by Balcan et al [1 7/. Those authors also present an algorithm 
with roughly the same query complexity for this problem. However, their algorithm is specifically designed for 
the linear threshold problem. Remarkably, near-optimal query complexity is achieved in polynomial-time by the 
general-purpose modified SGBS algorithm. 

^Note that the factor ^1 — — '^^ ^ in the exponential rate parameter A is a positive constant stiictly less than 1. For a noise 

level a this factor is maximized by a value /3 £ (a, 1/2) which tends to (1/2 + a)/2 as a tends to 1/2. 
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IV. Agnostic GBS 

So far we have assumed that the correct hypothesis h* is in H. In this section we drop this assumption and 
consider agnostic algorithms guaranteed to find the best hypothesis in H even if the correct hypothesis h* is not in 
H and/or the assumptions of Theorem [T] or [4] do not hold. The best hypothesis in H can be defined as the one that 
minimizes the error with respect to a probability measure on X, denoted by Px, which can be arbitrary. This notion 
of "best" commonly arises in machine problems where it is customary to measure the error or risk with respect 
to a distribution on X. A common approach to hypothesis selection is empirical risk minimization (ERM), which 
uses queries randomly drawn according to and then selects the hypothesis in 7i that minimizes the number of 
errors made on these queries. Given a budget of n queries, consider the following agnostic procedure. Divide the 
query budget into three equal portions. Use GBS (or NGBS or modified SGBS) with one portion, ERM (queries 
randomly distributed according Px) with another, and then allocate the third portion to queries from the subset 
of X where the hypothesis selected by GBS (or NGBS or modified SGBS) and the hypothesis selected by ERM 
disagree, with these queries randomly distributed according to the restriction of Px to this subset. Finally, select the 
hypothesis that makes the fewest mistakes on the third portion as the final choice. The sample complexity of this 
agnostic procedure is within a constant factor of that of the better of the two competing algorithms. For example, 
if the conditions of Theorems [T] or |4] hold, then the sample complexity of the agnostic algorithm is proportional 
to log N. In general, the sample complexity of the agnostic procedure is within a constant factor of that of ERM 
alone. We formalize this as follows. 

Lemma 4: Let Px denote a probability measure on X and for every h ^ % let R{h) denote its probability of 
error with respect to Px- Consider two hypotheses /ii, /12 G ^ cind let A d X denote the subset of queries for which 
hi and disagree; i.e., hi{x) ^ h2{x) for all a; G A. Suppose that m queries are drawn independently from Pa, 
the restriction of Px to the set A, let i?A(^i) cind -Ra(^2) denote average number of errors made by hi and /12 
on these queries, let Rj^{hi) =E[-Ra(^i)] and RA{h2) = E[_Ra(^2)], and select h = aTgimii{R^{hi), RA{h2)}. 
Then R(h) > xmn{R{hi) , R{h2)} with probability less than e~^MR^{hi)~R^{h2)\'' 

Proof Define 5 :— -Ra(^i) — ^^a(^2) and let 5 — RA{hi) — -Ra(/^2)- By Hoeffding's inequality we have 
5 e [5-e,6 + e] with probability at least 1 - 2e-"l'l'/2. it follows that P(sign(?) ^ sign((5)) < 2e-"l'^l'/2. For 
example, if 5 > then since ¥{6 < S - e) < 2e-'"''/2 we may take S = e to obtain P(? < 0) < 2e-"'''/2. The 
result follows since hi and /12 agree on the complement of A. ■ 
Note that there is a distinct advantage to drawing queries from Pa rather than Px, since the error exponent is 
proportional to |Pa(^i) — ^a(^2)P which is greater than or equal to \R{hi) ~ i?(ft,2)P- Now to illustrate the 
idea, consider an agnostic procedure based on modified SGBS and ERM. The following theorem is proved in the 
Appendix. 

Theorem 5: Let Px denote a measure on X and suppose we have a query budget of n. Let hi denote the 
hypothesis selected by modified SGBS using n/3 of the queries and let /12 denote the hypothesis selected by ERM 
from n/3 queries drawn independently from Px. Draw the remaining n/3 queries independently from P\, the 
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restriction of Px to the set A on which hi and /i2 disagree, and let R/^{hi) and R/^{h2) denote the average 
number of errors made by hi and /12 on these queries. Select h = arg inin{i?A(^i), -Ra(^2)}- Then, in general, 

E[i?(/i)] < min{E[i?(/ii)],E[i?(/i2)]} + VW^ : 

where R{h), h £ %, denotes the probability of error of h with respect to Px- Furthermore, if the assumptions of 
Theorem |4] hold and h* £ H, then 

V(h^h*) < Are^^"/3 + 2e-"|i-2"l'/6 , 

where a is the noise bound. 

Note that if the assumptions of Theorem |4] hold, then the agnostic procedure performs almost as well as modified 
SGBS alone. In particular, the number of queries required to ensure that ¥{h =^ h*) < 5 is proportional to log ^; 
optimal up to constant factors. Also observe that the probability bound implies the following bound in expectation: 

E[R(h)] < R{h*)+ Ne-^"/^ + 2e-"|i-2"lV6 

< R{h*) + Ne~^" , 

where C > is a constant depending on A and a. The exponential convergence of the expected risk is much faster 
than the usual parametric rate for ERM. 

If the conditions of Theorem [4] are not met, then modified SGBS (alone) may perform poorly since it might 
select inappropriate queries and could even terminate with an incorrect hypothesis. However, the expected error of 
the agnostic selection h is within \/3/n of the expected error of ERM, with no assumptions on the underlying 
distributions. Note that the expected error of ERM is proportional to n^^/^ in the worst-case situation. Therefore, 
the agnostic procedure is near-optimal in general. The agnostic procedure offers this safeguard on performance. 
The same approach could be used to derive agnostic procedures from any active learning scheme (i.e., learning 
from adaptively selected queries), including GBS or NGBS. We also note that the important element in the agnostic 
procedure is the selection of queries; the proposed selection of h based on those queries is convenient for proving 
the bounds, but not necessarily optimal. 

V. Applications OF GBS 

In this section we examine a variety of common situations in which the neighborliness condition can be verified. 
We will confine the discussion to GBS in the noise-free situation, and analogous results hold in the presence of 
noise. For a given pair {X,H), the effectiveness of GBS hinges on determining (or bounding) c* and establishing 
that (XjH) are neighborly. Recall the definition of the bound c* from ([T|| and that N — iH], the cardinality of T-L. 
A trivial bound for c* is 

max h{A) P{A) < 1 - , 
''^"^ AeA 

which is achieved by allocating N^^ mass to each hypothesis and evenly distributing mass to set(s) A where 
h{A) = +1 and on set(s) A where h{A) = —1. Non-trivial coherence bounds are those for which there exists 
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a P and < c < 1 that does not depend on N such that 

< c 



max 

hen 



Y^h{A)P{A) 

AeA 

The coherence parameter c* is analytically determined or bounded in several illustrative applications below. We 
also note that it may be known a priori that c* is bounded far way from 1. Suppose that for a certain P on X 
(or A) the absolute value of the first-moment of the correct hypothesis (w.r.t. P) is known to be upper bounded 
by a constant c < 1. Then all hypotheses that violate the bound can be eliminated from consideration. Thus the 
constant c is an upper bound on c* . Situations Uke this can arise, for example, in binary classification problems with 
side/prior knowledge that the marginal probabilities of the two classes are somewhat balanced. Then the moment 
of the correct hypothesis, with respect to the marginal probability distribution on X, is bounded far away from 1 
and —1. 

The neighborly condition can be numerically verified in a straightforward fashion. Enumerate the equivalence 
sets in as Ai, . . . ,Am- To check whether the /c-neighborhood graph is connected, form an M x M matrix R^. 
whose i,j entry is 1 if Ai and Aj are k-neighbors and otherwise. Normalize the rows of Rk so that each sums 
to 1 and denote the resulting stochastic matrix by Qk- The fc-neighborhood graph is connected if and only if there 
exists an integer £, < £ < M, such that contains no zero entries. This follows from the standard condition for 
state accessibility in Markov chains (for background on Markov chains one can refer to the textbook by Bremaud 
pT)). The smallest k for which the fc-neighborhood graph is connected can be determined using a binary search 
over the set 1, . . . , M, checking the condition above for each value of k in the search. This idea was suggested 
to me by Clayton Scott. Thus, the neighborly condition can be verified in polynomial-time in M. Alternatively, in 
many cases the neighborly condition can be verified analytically, as demonstrated in following applications. 



A. One Dimensional Problems 

First we show that GBS reduces to classic binary search. Let T-L — {ft-i, . . . , h^} be the collection of binary-valued 
functions on A" = [0, 1] of the following form, hi{x) := sign — jqip^j for i = 1, . . . ,N (and sign(O) := +1). 
Assume that h* e %. 

First consider the neighborly condition. Recall that A is the smallest partition of X into equivalence sets induced 
by %. In this case, each A is an interval of the form Ai = FTT-'' * = li • ■ • : Observe that only a single 

hypothesis, hi, has different responses to queries from Ai and Ai+i and so they are 1-neighbors, for i = 1, . . . , A^— 1. 
Moreover, the 1 -neighborhood graph is connected in this case, and so {X ,%) is 1-neighborly. 

Next consider coherence parameter c*. Take P to be two point masses at a; = and x = 1 of probability 1/2 
each. Then \Y.AeAKA)P{A)\ = for every h e U, since h{{)) = -1 and h{l) = +1. Thus, c* = 0. Since 
c* — Q and fc = 1, we have a = 2/3 and the query complexity of GBS is proportional to logA^ according to 
Theorem [T] The reduction factor of 2 /3, instead of 1 /2, arises because we allow the situation in which the number 
of hypotheses may be odd (e.g., given three hypotheses), the best query may eliminate just one). If A^ is even, then 
the query complexity is log2(Af), which is information-theoretically optimal. 
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Now let X = [0,1] and consider a finite collection of hypotheses Ti. = {hi{x)}fLi, where hi takes the value 
+1 when X E [a^, bi), for a pair < < bi < 1, and —1 otherwise. Assume that h* E %. The partition A again 
consists of intervals, and the neighborly condition is satisfied with fc = 1. To bound c* , note that the minimizing 
P must place some mass within and outside each interval [0^,6^). If the intervals all have length at least £ > 0, 
then taking P to be the uniform measure on [0, 1] yields that c* < 1 — 2£, regardless of the number of interval 
hypotheses under consideration. Therefore, in this setting Theorem [T] guarantees that GBS determines the correct 
hypothesis using at most a constant times log N steps. 

However, consider the special case in which the intervals are disjoint. Then it is not hard to see that the best 
allocation of mass is to place mass in each subinterval, resulting in c* = 1 — 2N^^. And so. Theorem [l]only 
guarantees that GBS is will terminate in at most N steps (the number of steps required by exhaustive linear search). 
In fact, it is easy to see that no procedure can do better than linear search in this case and the query complexity of 
any method is proportional to N . However, note that if queries of a different form were allowed, then much better 
performance is possible. For example, if queries in the form of dyadic subinterval tests were allowed (e.g., tests 
that indicate whether or not the correct hypothesis is +l-valued anywhere within a dyadic subinterval of choice), 
then the correct hypothesis can be identified through [log2 queries (essentially a binary encoding of the correct 
hypothesis). This underscores the importance of the geometrical relationship between X and % embodied in the 
neighborly condition and the incoherence parameter c* . Optimizing the query space to the structure of Ti. is related 
to the notion of arbitrary queries examined in the work of KuUcarni et al |25|, and somewhat to the theory of 



compressed sensing developed by Candes et al |39| and Donoho pO). 



B. Multidimensional Problems 

Let % — be a collection of multidimensional threshold functions of the following form. The threshold 

of each hi determined by (possibly nonlinear) decision surface in rf-dimensional Euclidean space and the queries 
are points in X It suffices to consider linear decision surfaces of the form 

K{x) ■.^?.\gn{{a^,x) +bi), (5) 

where E M'^, ||ai||2 = 1, the offset 6^ G M satisfies \bi\ < b for some constant b < 00, and {ai,x) denotes the 
inner product in M"^. Each hypothesis is associated with a halfspace of M''. Note that hypotheses of this form can 
be used to represent nonlinear decision surfaces, by first applying a mapping to an input space and then forming 
linear decision surfaces in the induced query space. The problem of learning multidimensional threshold functions 



arises commonly in computer vision (see the review of Swain and Strieker |12| and applications by Geman and 
Jedynak p3| and Arkin et al fTTi), image processing studied by Korostelev and Kim pO| , pT| , and active learning 
research; for example the investigations by Freund et al fTSl, Dasgupta p6) , Balcan et al |17|, and Castro and 
Nowak imi. 



First we show that the pair {M.'^,H) is 1-neighborly. Each A E A is a polytope in M.'^. These polytopes are 
generated by intersections of the halfspaces corresponding to the hypotheses. Any two polytopes that share a 
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common face are 1-neighbors (the hypothesis whose decision boundary defines the face, and its complement if 
it exists, are the only ones that predict different values on these two sets). Since the polytopes tessellate M'', the 
1-neighborhood graph of A is connected. 

We next show that the the coherence parameter c* = 0. Since the offsets of the hypotheses are all less than 
b in magnitude, it follows that the distance from the origin to the nearest point of the decision surface of every 
hypothesis is at most b. Let Pr denote the uniform probability distribution on a ball of radius r centered at the 
origin in M''. Then for every h of the form (jsj) there exists a constant C > (depending on b) such that 

C 



Y.KA)PM) 



h{x) dPr{x) 



< 



and limr^oolJ-gid h{x) dPr{x)\ = 0. Therefore c* = 0, and it follows from Theorem [T| guarantees that GBS 
determines the correct multidimensional threshold in at most [ ;og°|/2) l'^g2 ^1 steps. To the best of our knowledge 
this is a new result in the theory of learning multidimensional threshold functions, although similar query complexity 
bounds have been established for the subclass of linear threshold functions with 6 = (threshold boundaries passing 
through the origin); see for example the work of Balcan et al |i7|. These results are based on somewhat different 
learning algorithms, assumptions and analysis techniques. 

Observe that if "H is an e-dense (with respect Lesbegue measure over a compact set in M'*) subset of the continuous 
class of threshold functions of the form then the size of the H satisfies log {Hi = log oc d\oge~^. Therefore 
the query complexity of GBS is proportional to the metric entropy of the continuous class, and it follows from 
the results of Kulkarni et al |25 1 that no learning algorithm exists with a lower query complexity (up to constant 
factors). Furthermore, note that the computational complexity of GBS for hypotheses of the form (|5]l is proportional 
to the cardinality of A, which is equal to the number of polytopes generated by intersections of half-spaces. It is a 
weU known fact (see Buck |42|) that \A\ = J^Uo i'j) °^ Therefore, GBS is a polynomial-time algorithm for 
this problem. In general, the cardinality of A could be as large as 2^. 

Next let Ti. again be the hypotheses of the form (jsjl, but let X [—1, 1]'', instead of all of M.'^. This constraint 
on the query space affects the bound on the coherence c*. To bound c*, let P be point masses of probability 2^'* 
at each of the 2"^ vertices of the cube [—1,1]'' (the natural generalization of the P chosen in the case of classic 



binary search in Section V-A above). Then jX^Ae^ '^l^) -P(^) | ~ \fx^(^^'^-^(^)\ — ^ ^ 2^''+^ for every 
h E H, since for each h there is at least one vertex on where it predicts +1 and one where it predicts —1. Thus, 
c* < 1 — 2^'*+^. We conclude that the GBS determines the correct hypothesis in proportional to 2'^logiV steps. 
The dependence on 2'' is unavoidable, since it may be that each threshold function takes that value +1 only at one 
of the 2'* vertices and so each vertex must be queried. A noteworthy special case is arises when b = (i.e., the 
threshold boundaries pass through the origin). In this case, with P as specified above, c* = 0, since each hypothesis 
responds with +1 at half of the vertices and —1 on the other half. Therefore, the query complexity of GBS is at 
most \ jog°|/2) log2 ^1' independent of the dimension. As discussed above, similar results for this special case have 



been previously reported based on different algorithms and analyses; see the results in the work of Balcan et al 1 17 1 
and the references therein. Note that even if the threshold boundaries do not pass through the origin, and therefore 
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the number of queries needed is proportional to logA^ so long as |6| < 1. The dependence on dimension d can 
also be eliminated if it is known that for a certain distribution P on X the absolute value of the moment of the 
correct hypothesis w.rt. P is known to be upper bounded by a constant c < 1, as discussed at the beginning of this 
section. 

Finally, we also mention hypotheses associated with axis-aligned rectangles in [0,1]'^, the multidimensional 
version of the interval hypotheses considered above. An axis-aligned rectangle is defined by its boundary coordinates 
in each dimension, {aj, bj}j^^, < aj < bj < 1. The hypothesis associated with such a rectangle takes the value 
+ 1 on the set {a; e [0, 1]"^ : aj < Xj < bj, j = 1, . . . ,d} and —1 otherwise. The complementary hypothesis 
may also be included. Consider a finite collection H of hypotheses of this form. If the rectangles associated with 
each h E H have volume at least ly, then by taking P to be the uniform measure on [0, 1]'' it follows that the 
coherence parameter c* < I — 2^ for this problem. The cells of partition A of [0, 1] associated with a collection of 
such hypotheses are rectangles themselves. If the boundaries of the rectangles associated with the hypotheses are 
distinct, then the 1 -neighborhood graph of A is connected. Theorem [T] implies that the number of queries needed 
by GBS to determine the correct rectangle is proportional to logiV/log((l — i^)^^)- 

C. Discrete Query Spaces 

In many situations both the hypothesis and query spaces may be discrete. A machine learning application, for 
example, may have access to a large (but finite) pool of unlabeled examples, any of which may be queried for a label. 
Because obtaining labels can be costly, "active" learning algorithms select only those examples that are predicted 
to be highly informative for labeling. Theorem [T] applies equally well to continuous or discrete query spaces. For 
example, consider the linear separator case, but instead of the query space R'^ suppose that A" is a finite subset of 
points in M'^. The hypotheses again induce a partition of X into subsets A{X, H), but the number of subsets in the 
partition may be less than the number in A{M.'^,'H). Consequently, the neighborhood graph of A{X,'H) depends 
on the specific points that are included in X and may or may not be connected. As discussed at the beginning of 
this section, the neighborly condition can be verified in polynomial-time (polynomial in \A\ < \X\). 

Consider two illustrative examples. Let "H be a collection of linear separators as in (|5]l above and first reconsider 
the partition A{M.'',T-L). Recall that each set in A{M.'^,H) is a polytope. Suppose that a discrete set X contains at 
least one point inside each of the polytopes in A{M.'^,7i). Then it follows from the results above that {X,H) is 1- 
neighborly. Second, consider a simple case in d = 2 dimensions. Suppose X consists of just three non-colinear points 
{xi, X2, x^} and suppose that H is comprised of six classifiers, {hf , , /i^, h'^ , h^}^ satisfying hf{xi) = +1, 
hf{xj) = — 1, j 7^ i , i = 1, 2, 3, and = —hf, i = 1, 2, 3. In this case, AiX^V.) — {{a^i}, {2^2}, {xs}} and the 
responses to any pair of queries differ for four of the six hypotheses. Thus, the 4-neighborhood graph of A{X, Ti) 
is connected, but the 1 -neighborhood is not. 

Also note that a finite query space naturally limits the number of hypotheses that need be considered. Consider 
an uncountable collection of hypotheses. The number of unique labeling assignments generated by these hypotheses 
can be bounded in terms of the VC dimension of the class; see the book by Vapnik for more information on VC 
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theory |43 1. As a result, it suffices to consider a finite subset of the hypotheses consisting of just one representative 
of each unique labeling assignment. Furthermore, the computational complexity of GBS is proportional to N \X\ 
in this case. 



VI. Related Work 

Generalized binary search can be viewed as a generalization of classic binary search, Shannon-Fano coding as 
noted by Goodman and Smyth |[T], and channel coding with noiseless feedback as studied by Horstein ||2]. Problems 
of this nature arise in many applications, including channel coding (e.g., the work of Horstein Q and Zigangirov 
Q), experimental design (e.g., as studied by Renyi ||6)), disease diagnosis (e.g., see the work of Loveland |7|), 
fault-tolerant computing (e.g., the work of Feige et al ||8]), the scheduling problem considered by Kosaraju et al 
computer vision problems investigated by Geman and Jedynak p3] and Arkin et al |14|), image processing 
problems studied by Korostelev and Kim pO| , |[TT], and active learning research; for example the investigations by 
Freund et al ||T5), Dasgupta |[16), Balcan et al {IT), and Castro and Nowak (31]. 

Past work has provided a partial characterization of this problem. If the responses to queries are noiseless, then 
selecting the sequence of queries from X is equivalent to determining a binary decision tree, where a sequence of 
queries defines a path from the root of the tree (corresponding to H) to a leaf (corresponding to a single element of 
H). In general the determination of the optimal (worst- or average-case) tree is NP-complete as shown by Hyafil and 
Rivest [[19|. However, there exists a greedy procedure that yields query sequences that are within a factor of log 
of the optimal search tree depth; this result has been discovered independently by several researchers including 
Loveland |j7), Garey and Graham |]7), Arkin et al pO| , and Dasgupta p6) . The greedy procedure is referred to here 
as Generalized Binary Search (GBS) or the splitting algorithm, and it reduces to classic binary search, as discussed 
in Section IV-AI 

The number of queries an algorithm requires to determine h* is called the query complexity of the algorithm. 
Since the hypotheses are assumed to be distinct, it is clear that the query complexity of GBS is at most N (because 
it is always possible to find query that eliminates at least one hypothesis at each step). In fact, there are simple 
examples (see Section |V-A| i demonstrating that this is the best one can hope to do in general. However, it is also 
true that in many cases the performance of GBS can be much better, requiring as few as log2(A^) queries. In 
classic binary search, for example, half of the hypotheses are eliminated at each step (e.g., refer to the textbook by 
Gormen et al pT[). Renyi first considered a form of binary search with noise |j5] and explored its connections with 
information theory ||6). In particular, the problem of sequential transmission over a binary symmetric channel with 
noiseless feedback, as formulated by Horstein |j2) and studied by Burnashev and Zigangirov (3) and more recently 
by Pelc et al p2) , is equivalent to a noisy binary search problem. 

There is a large literature on learning from queries; see the review articles by Angluin ||23|, p4) . This paper 
focuses exclusively on membership queries (i.e., an x E X is the query and the response is h*{x)), although 
other types of queries (equivalence, subset, superset, disjointness, and exhaustiveness) are possible as discussed by 
Angluin p3| . Arbitrary queries have also been investigated, in which the query is a subset of H and the output 
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is +1 if h* belongs to the subset and —1 otherwise. A finite collection of hypotheses H can be successively 
halved using arbitrary queries, and so it is possible to determine h* with log2 N arbitrary queries, the information- 
theoretically optimal query complexity discussed by KuUcarni et al |25|. Membership queries are the most natural 
in function learning problems, and because this paper deals only with this type we will simply refer to them as 
queries throughout the rest of the paper. The number of queries required to determine a binary-valued function in a 
finite collection of hypotheses can be bounded (above and below) in terms of a combinatorial parameter of {X, Ti) 
due to Hegediis | p6| (see the work of Hellerstein et al |27| for related work). Due to its combinatorial nature, 
computing such bounds are generally NP-hard. In contrast, the geometric relationship between X and % developed 
in this paper leads to an upper bound on the query complexity that can be determined analytically or computed in 
polynomial time in many cases of interest. 

The term GBS is used in this paper to emphasize connections and similarities with classic binary search, which is 
a special case the general problem considered here. Classic binary search is equivalent to learning a one-dimensional 
binary-valued threshold function by selecting point evaluations of the function according to a bisection procedure. 
Consider the threshold function ht{x) :— sign(a; — <) on the interval X := [0, 1] for some threshold value t E (0, 1). 
Throughout the paper we adopt the convention that sign(O) = +1. Suppose that t belongs to the discrete set 
iw+i' ■ ■ ■ ' TTTi} ^ denote the collection of threshold functions hi, . . . ,h]y. The value of t can then 

determined from a constant times logA^ queries using a bisection procedure analogous to the game of twenty 
questions. In fact, this is precisely what GBS performs in this case (i.e., GBS reduces to classic binary search 
in this setting). If iV = 2'" for some integer m, then each point evaluation provides one bit in the m-bit binary 
expansion of t. Thus, classic binary search is information-theoretically optimal; see the book by Traub, Wasilkowski 
and Wozniakowski p8) for a nice treatment of classic bisection and binary search. The main results of this paper 
generalize the salient aspects of classic binary search to a much broader class of problems. In many (if not most) 
applications it is unrealistic to assume that the responses to queries are without error. A form of binary search 
with noise appears to have been first posed by Renyi ||5j. The noisy binary search problem arises in sequential 
transmission over a binary symmetric channel with noiseless feedback studied by Horstein ||2j and Zigangirov Q, 
p9) . The survey paper by Pelc et al p2) discusses the connections between search and coding problems. In channel 
coding with feedback, each threshold corresponds to a unique binary codeword (the binary expansion of t). Thus, 
channel coding with noiseless feedback is equivalent to the problem of learning a one-dimensional threshold function 
in binary noise, as noted by Burnashev and Zigangirov ||3). The near-optimal solutions to the noisy binary search 
problem first appear in these two contexts. Discrete versions of Horstein's probabilistic bisection procedure Q were 
shown to be information-theoretically optimal (optimal decay of the error probability) in the works of Zigangirov 
and Burnashev Q, p9] . More recently, the same procedure was independently proposed and analyzed in the 
context of noise-tolerant versions of the classic binary search problem by Karp and Kleinberg pO| , which was 
motivated by applications ranging from investment planning to admission control in queueing networks. Closely 
related approaches are considered in the work of Feige et al |8|. The noisy binary search problem has found 
important applications in the minimax theory of sequential, adaptive sampling procedures proposed by Korostelov 
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and Kim 1 10 1, fTT) for image recovery and binary classification problems studied by Castro and Nowak | [3Tj . We 



also mention the works of Rivest et al | 32| , Spencer |33J and Aslam and Dhagat p5| , and Dhagat et al p4) , which 
consider adversarial situations in which the total number of erroneous oracle responses is fixed in advance. 

One straightforward approach to noisy GBS is to follow the GBS algorithm, but to repeat the query at each step 
multiple times in order to decide whether the response is more probably +1 or —1. This simple approach has been 
studied in the context of noisy versions of classic binary search and shown to perform significantly worse than 
other approaches in the work of Karp and Kleinberg pOj ; perhaps not surprising since this is essentially a simple 
repetition code approach to communicating over a noisy channel. A near-optimal noise-tolerant version of GBS 
was developed in this paper. The algorithm can be viewed as a non-trivial generalization of Horstein's probabilistic 
bisection procedure. Horstein's method relies on the special structure of classic binary search, namely that the 
hypotheses and queries can be naturally ordered together in the unit interval. Horstein's method is a sequential 
Bayesian procedure. It begins with uniform distribution over the set of hypotheses. At each step, it queries at the 
point that bisects the probability mass of the current distribution over hypotheses, and then updates the distribution 
according to Bayes rule. Horstein's procedure isn't directly applicable to situations in which the hypotheses and 
queries cannot be ordered togetherl, but the geometric condition developed in this paper provides similar structure 
that is exploited here to devise a generalized probabilistic bisection procedure. The key elements of the procedure 
and the analysis of its convergence are fundamentally different from those in the classic binary search work of 
Burnashev and Zigangirov |[3| and Karp and Kleinberg | |30| . 

VII. Conclusions and Possible Extensions 

This paper investigated a generalization of classic binary search, called GBS, that extends it to arbitrary query 
and hypothesis spaces. While the GBS algorithm is well-known, past work has only partially characterized its 
capabilities. This paper developed new conditions under which GBS (and a noise-tolerant variant) achieve the 
information-theoretically optimal query complexity. The new conditions are based on a novel geometric relation 
between the query and hypothesis spaces, which is verifiable analytically and/or computationally in many cases of 
practical interest. The main results are applied to learning multidimensional threshold functions, a problem arising 
routinely in image processing and machine learning. 

Let us briefly consider some possible extensions and open problems. First recall that in noisy situations it 
is assumed that the binary noise probability has a known upper bound a < 1/2. It is possible to accommodate 
situations in which the bound is unknown a priori. This can be accomplished using an NGBS algorithm in which the 
number of repetitions of each query, R, is determined adaptively to adjust to the unknown noise level. This procedure 
was developed by the author in p8) , and is based on a straightforward, iterated application of Chernoff's bound. 
Similar strategies have been suggested as a general approach for devising noise-tolerant learning algorithms (36). 
Using an adaptive procedure for adjusting the number of repetitions of each query yields an NGBS algorithm with 
query complexity bound proportional to log log the same order as that of the NGBS algorithm discussed 

above which assumed a known bound a. Whether or not the additional logarithmic factor can be removed if the 
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noise bound a is unknown is an open question. Adversarial noise models in which total number of errors is fixed 
in advance, like those considered by Rivest et al |32] and Spencer |33|, are also of interest in classic binary search 
problems. Repeating each query multiple times and taking the majority vote of the responses, as in the NGBS 
algorithm, is a standard approach to adversarial noise. Thus, NGBS provides an algorithm for generalized binary 
search with adversarial noise. 

Finally, we suggest that the salient features of the GBS algorithms could be extended to handle continuous, 
uncountable classes of hypotheses. For example, consider the continuous class of halfspace threshold functions 
on K''. This class is indexed parametrically and it is possible to associate a volume measure with the class (and 
subsets of it) by introducing a measure over the parameter space. At each step of a GBS-style algorithm, all 
inconsistent hypotheses are eliminated and the next query is selected to split the volume of the parameter space 
corresponding to the remaining hypotheses, mimicking the splitting criterion of the GBS algorithm presented here. 

Acknowledgements. The author thanks R. Castro, A. Gupta, C. Scott, A. Singh and the anonymous 
reviewers for helpful feedback and suggestions. 

Vin. Appendix 

A. Proof of Lemma |2] 

First we derive the precise form of pi,p2, ... is derived as follows. Let 6i = {i + J2hPii^) ^i(^))/2, the weighted 
proportion of hypotheses that agree with yi. The factor that normaUzes the updated distribution in (|2| is related to Si 
as follows. Note that /jd-' - = E„:..(„)=-i + E/.:..(^.)=i " /?) = 

(1- (5,0/? + 5.(1-/3). Thus, 

^(l-z>('i))/2(l _ _0)(l+z,(/i))/2 

^■'-^^^'^ ^ ^'^('^ (l-<5.)/3 + 5.(1-/3) 
Denote the reciprocal of the update factor for pi+i{h*) by 

(1-50/3 + 5,(1-/3) 

where Zi{h*) = h*{xi)yi, and observe that pi^i{h*) = pi{h*)/^i. Thus, 

Q+i ^ {l-p,{h*)hi)p,{h*) ^ -U-P^{h*) 
k(/i*)/7,(1-p,(/i*)) l-p^{h*) 

We prove that E[Ci+i|pi] < Ci by showing that E[7i|pi] < 1. To accomplish this, we will let pi be arbitrary. 

For every A £ A and every ft. G "H let h{A) denote the value of h on the set A. Define 5^ = 

(1 + J2hPiW^i^)) f^' '^he proportion of hypotheses that take the value +1 on A. Let Ai denote that set that 
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Xi is selected from, and consider the four possible situations: 





= +1, Vi 


= +1 : 


7» = 


(1- 


-5+ 


-/3) 




1-/3 






= +1, Vi 


= -1 : 


l^ = 


^1 


P+{1-S\,)(l- 


-/3) 




P 




h*{xi) 


= -1, Vi 


= +1 : 


li = 






-/3) 




p 




h*{x,) 


= -1, y^ 


= -1 : 


li = 




/3+(l-5+^)(l- 


-/3) 




1-/3 





To bound E[7i|pj] it is helpful to condition on Ai. Define q 



^'y^\p^,A,] = 



(1-*) H ■ 



1-/3 



^ h*{xi)). If /i*(AO = +1, then 



= '5+ +(1-5+) 

Define 7+(A,):=<5X + (l-'5X)'^'^ 

E[7,;b„A,] = (1-5+) + ,5+ 



li{l-q,) , 9.(1-/3) 



1-/3 ^ /3 



1-/3 /3 

Similai-ly, if /i*(Aj) = -1, then 
/?(l-g,) , q,{l-P) 



1-/3 



By assumption q,<a<l/2, and since a < /3 < 1/2 the factor ^^jEf^ + 2^^^ < + < 1 (strictly 

less than 1 if /3 > a). Define 

_ /3(l-a) a(l-/3) 

to obtain the bounds 



lt{A^) < 5+ +(l-<5+)(l-eo) 
J- (A,) < 5+ (l-eo) + (l-^X) 



(7) 
(8) 



Observe that for every A we have < 5^ < 1, since at least one hypothesis takes the value —1 on A and p{h) > 
for all h £ v.. Therefore both ■y^{Ai) and ■y^ (Ai) are less or equal to 1, and it follows that E[7i|pi] < 1 (and 
strictly less than 1 if /3 > a). ■ 

B. Proof of Theorem [4] 

The proof amounts to obtaining upper bounds for yf{Ai) and 7" (A;), defined above in (j7| and (jsji. Consider two 
distinct situations. Define hi :— min^g_4 |VF(pi, A)|. First suppose that there do not exist neighboring sets A and 
A' with W{pi, A) > bi and W{pi,A') < —bi. Then by Lemma[T| this implies that bi < c*, and according the query 
selection step of the modified SGBS algorithm, Ai — argmin^ |M^(pi, A)|. Note that because |Ty(pi, Ai)| < c*, 
(l-c*)/2 < S^^ < (1 + c*)/2. Hence, both 7+(A,) and 7~(A,) are bounded above by 1 - £0(1 - c*)/2. 

Now suppose that there exist neighboring sets A and A' with W{pi,A) > bi and W{pi,A') < —bi. Recall 
that in this case Ai is randomly chosen to be A or A' with equal probability. Note that > {I + bi)/2 and 
5% < (1 - bi)/2. If h*{A) = h*{A') = +1, then applying (j?} results in 

E[y,\p,,Ae{A,A'}] < i(i + l_^ + l±^(i_eo)) = i(2-eo^) < l-£o/4, 



;{2-eoil + S+-5+)) 



= l-^{l + S+-6+) < l-eo/2 
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since b, > 0. Similai-ly, if h*{A) = h*{A') = -1, then ^ yields E[7,|p,, e {A,A'}] < l-£o/4. lfh*{A) = -1 
on A and h*{A) = +1, then applying ([8]l on A and (j?]) on A' yields 

E[^,\p,,A&{A,A}] < ^{5+il-eo) + {l-5^) + 6+, + il-S+){l-eo)) 

- \{l-S+ + S++{l~eo){l + S+-S+)) 
1 

2 

since < (J;^ - (5^, < 1. The final possibility is that h*{A) = +1 and h*{A') = -1. Apply (j7| on A and ([sjl on 
A to obtain 

E[7,|p„ A, G A'}] < ^ (<5+ + (1 - <5+)(l - eo) + <5+ (1 - So) + (1 - <5+ )) 

- l{l + S+-S++{l-eo){l-6++6+)) 

Next, use the fact that because A and A' are neighbors, 6^ — 6^, — Pi{h*) — pi{—h*); if ~h* does not belong to 
n, then pr{-h*) = 0. Hence, 

= ^{l+p,{h*)-p,{~h*) + (1 - eo)(l 

< i(l+p,(r) + (l-eo)(l-p.(/i*))) = , 

since the bound is maximized when pi{—h*) = 0. Now bound E[7,;|pi] by the maximum of the conditional bounds 
above to obtain 

E[7.N < max{l-^(l~K(/i*)),l-|,l-(l-c*)|} , 
and thus it is easy to see that 



E 



Ik 



C. Proof of Theorem [5] 



E[7,|p,]-K(fe*) ^ 
l-p^{h*) 



l-min{|(l-c*),f} 



First consider the bound on W[R{h)]. Let (5„ = 2 e^"'^^*^''!^"^^^''^^!^/^ and consider the conditional expectation 
E[i?(/i)|/ii, /12]; i.e., expectation with respect to the n/3 queries drawn from the region A, conditioned on the 2n/3 
queries used to select hi and ft,2. By Lemma |4] 

E[i?(/i)|/ii,/i2] < (l-(5„) min{i?(/ii),i?(/i2)} + (Jn max{i?(/ii),i?(/i2)} , 

= min{i?(/ii),i?(/i2)} + (5„ [max{i?(/ii),i?(ft2)}-min{i?(/ii),i?(/i2)}] , 

= min{i?(/ii),i?(;i2)} + 5„|i?(/ii)-i?(/i2)| , 

= m:in{R{hi),R[h2)} + 2 - i?(/i2)| e-"l^^(''i)-^^(''^)l'/^ , 

< m\n{R{hi),R{h2)} + 2|i?(/ii) - i?(/i2)| e-'^l^e^^^-^C'^'l'/^ , 
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where the last inequaUty follows from the fact that \R{hi) — R{h2)\ < |i?A(^i) ^ -Ra(^2)|- 

^ 2 / Q ^ — 

The function 2ue^^" attains its maximum at u = w ^ < yl/C' ™d therefore 

E[R(h)\hi,h2] < mm{R{hi),R{h2)} + y^S/n . 

Now taking the expectation with respect to hi and /i2 (i.e., with respect to the queries used for the selection of hi 
and ft,2) 

E[R(h)] < E[nim{R{hi),R{h2)}] + ^/3/^ , 

< min{E[R{hi)],E[R{h2)]} + ^/sjn , 

by Jensen's inequality. 

Next consider the bound on P{hi ^ h*). This also follows from an application of Lemma |4] Note that if the 
conditions of Theorem |4]hold, then P(/ii ^ h*) < iVe-^"/^ Purthemore, if hi = h* and ft-a ^ h* , then 
I^a(^i) — -Ra(^2)| > |1 — 2q;|. The bound on V{h ^ h*) follows by applying the union bound to the events 
hi = h* and R(h) > mm{R{hi), R{h2)}. ■ 
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